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Abstract. A liquid foam in contact with a solid surface forms a two-dimensional foam on the surface. We 
derive the equilibrium equations for this 2D foam when the solid surface is curved and smooth, generalising 
the standard case of flat Hele Shaw cells. The equilibrium conditions at the vertices in 2D, at the edges in 
3D, are invariant by conformal transformations. Regarding the films, conformal invariance only holds with 
restrictions, which we explicit for 3D and flat 2D foams. Considering foams confined in thin interstices 
between two non parallel plates, normal incidence and Laplace's law lead to an approximate equation 
relating the plate profile to the conformal map. Solutions are given for the logarithm and power laws in 
the case of constant pressure. The paper concludes on a comparison with available experimental data. 
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1 Introduction 

Many applications of liquid foams, from dish washing to 
shaving cream, container decontamination, penetration and 
dispersion through porous media like rocks, etc., involve 
contacts with solid surfaces. 

One of the most common settings, for both experi- 
mental and theoretical investigations, is the Hele Shaw 
cell consisting of two flat and parallel plates separated 
by a small distance. The foam, confined in the interstice, 
may often be considered as quasi bi-dimensional; the equi- 
librium equations are 2D versions of Plateau's equations 
|ll2j . Even in this simple setting, the predictions of straight- 
forward dimensional reduction may fail when the thickness 
of the bubbles are not much smaller than their diameter 
in the plane of the plates. Then instabilities occur |i3i4i5) . 

In the present paper, we consider standard dry liq- 
uid foams in contact with smooth curved solid surfaces. 
Later, we will come back to settings similar to the Hele 
Shaw cell. For the solid surface — the container wall in 
typical applications — , the main assumptions are rigidity 
and smoothness. Contact with non rigid walls such as elas- 
tic membranes or an interface with another fluid medium 
must be treated differently; in a number of cases, at least 
for statics, flexible surfaces may be included as part of 
the foam with a surface tension differing from the rest of 
the foam. Smoothness rules out corners, wedges as well as 
dirty or fractal surfaces. Usually, the films are pinned by 
such irregularities. 

^ CNRS UMR 8089 



On a smooth and clean surface S, the contacts of the 
real 3D foam form a 2D foam in S. Upon meeting the 
boundary surface, the thin films, or cell interfaces, define 
edges in S, the Plateau borders (edges in 3D) end up onto 
vertices in S and the cells are incident to faces, the 2D 
cells in S. The films are free to move along the surface, 
with viscous dissipation in dynamic conditions and full 
relaxation to equilibrium in static conditions. 

In recent experiments, the simple parallelepiped geom- 
etry was modified. Whereas the focus for hemispherical, 
yet parallel, plates was on topology Drenckhan et al. 
[7] tried several non parallel chambers and noticed that 
each resulting 2D foam could be related to the reference 
regular honeycomb pattern by a conformal map. The dom- 
inant constraint was volume conservation. A question ad- 
dressed in the present paper is whether the shape of the 
films, that is, the local geometry of the foam, can also be 
related to the conformal map and how. 

Conformality means that the transformation locally re- 
duces to a simple isotropic dilatation, preserving angles. 
However, in 2D, general conformal transformations are not 
compatible with Laplace law (at least, the 2D version). 
The only compatible conformal transformations (mapping 
circle arcs to circle arcs) are the Moebius transformations 
(homographies) |8I9I1(J| . In the case of algebraic functions 
(with exponents other than ±1), or the logarithm, or any 
other holomorphic function, is conformality only a long 
wavelength approximation ? or is there a correction to 
the 2D equilibrium ? 

On one hand, the equilibrium conditions at the Plateau 
borders — of dimension 2 less than the embedding space. 
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the edges in 3D, the vertices in 2D — turn out to be in- 
variant by conformal mappings. This result is exact within 
the standard models of dry foams in any dimension. On 
the other hand, regarding the films — bubbles interfaces, 
of co-dimension 1 — in thin cells, we will show that incor- 
porating the curvature in the third direction restores com- 
patibility with Laplace, implying a relation between the 
profile and the map. The case of constant pressure will be 
solved explicitly. Comparison with constant volume set- 
tings, including experimental ones, will also be done. 

Outline. After recalling some basics on foam equilib- 
rium and contact with a solid boundary in Section |21 
we derive, in Section |31 the equilibrium equations at the 
boundary in the case of normal incidence. Section^lestab- 
lishes conformal invariance for the co-dimension 2 Plateau 
borders equilibrium. Section treats the question of con- 
formal mappings for the films and solves specific examples 
to lowest order in non parallelism. Section El concludes. 
Appendix IXIandlBlrecall elements of geometry for non ex- 
perts. |0 contains some extensions dealing with the cases 
where the surface tension along the solid varies from cell 
to cell, so that the contact angle is not right. 



2 Foam equilibrium 

In equilibrated dry foams, the bubbles are separated by 
films obeying Laplace- Young's law: 

AP + 2jH = (1) 

where AP = P2 — Pi is the pressure drop, 7 the surface 
tension and H the mean curvature of the film. In static 
conditions and in absence of applied field, the pressure in- 
side each bubble is constant so that the films are constant 
mean curvature surfaces. 

The films and cells meet three by three at edges in a 
manner satisfying Plateau's laws : 

E ■=! 7,b, = (2) 
E -=1 IsH, = (3) 

The first law asserts that the sum of the forces on any 
edge is zero; indeed the force per unit edge length exerted 
by the film J^j is the surface tension times the unit vec- 
tor hj tangent to the film and normal to the junction. In 
general situations, the surface tension may have different 
values 7j on different sides. 

The second law Q is the sum of the pressure drops 
along a small circular path encircling the edge. By (O, 
this sum of differences is zero as long as the path is closed: 
(Pi - P2) + {P2 - P3) + {P3 Pi) = 0- Notice the sign 
convention implicitly assumed under Q: given a sense of 
circulation, the pressure difference in going from bubble 1 
to bubble 2 is defined as P2 — Pi (head minus tail) . By , 
this locally fixes the orientation of the three films around 
any given edge or Plateau border. The same film J- may 
appear with different signs in sums around different edges. 



2.1 Contact with a boundary 

The foam is in contact with a solid wall, not necessarily 
flat but smooth and clean. The films or interfaces can 
freely slide along the surface so as to relax to equilibrium. 

With these assumptions, an interface i has a contact 
with the wall only in non wetting conditions; then the 
contact angle 9 satisfies Young's law: 

7-, cos 6 + jsi- 1S2 = (4) 

where 1 and 2 label the two fluids (gas in normal foams) 
on both sides of the film. 

In the case of a film — a double locally symmetric 
interface — the two bubbles contain the same gas so that 
7is = 72s and cos 6* = 0. This condition of normal in- 
cidence is typical of soap froth; it is also valid when the 
soap solution wets the solid surface, provided the wetting 
films on both sides have similar surface tension. 

In the standard Hele-Shaw experimental setup, when 
the glass plates separation h is small, smaller than any 
typical length (curvature radius or cell diameter) of the 
2D foam, the curvature k± of the films in the direction 
normal to the plates is approximately zero. Then, the 2D 
foam satisfies equilibrium equations similar to (|2J), Q with 
quantities adapted to 2D: line tension 7^, edge tangent 
vector hj and curvature Hj — (2rj)~^, rj being the signed 
curvature radius of film i. In flat equilibrated 2D foams, 
the edges are circle arcs. 

Our purpose, next, is to consider non planar glass sur- 
faces. 

3 Equilibrium conditions at a boundary 

We derive equilibrium equations, analogous to Plateau's 
laws, valid for the 2D foam formed by the contacts of the 
films with the solid boundary. 

A Plateau border is a line -0 at the junction of three 
films J-j,i — 1,2,3. We call this a triple. These films 
and border are in contact with a solid smooth surface S. 
The contact of each film J-j gives rise to a line ^pj in S] 
these three lines emanate from a vertex p of the 2D foam 
that is the contact of the border -;/; with S (fig. Un- 
less otherwise specified, all the curves are parametrised 
by arc length s with s = at p. Then, for any curve Lp, 
T — ip = j^(p is the unit vector field tangent to (p. 

Assume normal incidence on S as implied by equilib- 
rium for real soap froth on clean surfaces. Then At p, 
common to Ti, J-2, and S , the following equations are 
verified: 

Ei=i 7jTj = 0, (5) 
Ei=i7ifcs(¥'i,'5) = 0, (6) 

where kg{ipj,S) is the geodesic curvature of ipj inS. 

The geodesic curvature is defined in several books (e.g. 
|11I12| ) and in Appendix A. 

When all the film tensions are equal, 7^ — 7, the tri- 
angles representing the vector sums in |j2J) and jSJ) are 
equilateral and the Plateau angles are 27r/3. 

For the proofs, we use two lemmas. 
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Fig. 1. The neighbourhood of a vertex p of a 2D foam: three 
films and a Plateau border meet at p on the solid surface S. 

3.1 Two lemmas 

The following lemmas are valid in full generality. The nor- 
mal incidence assumption will only be used to prove the 
main equations Q, I©. 

Lemma 1 Let J- ^ S he two smooth embedded surfaces in- 
tersecting transversely^ on a line (p. Then the geodesic and 
normal curvatures satisfy the following relation at any in- 
tersection point x: 

( ^3(95,5) \ _ ( COS0 -sin6'\ / kg{ip,T)\ , . 

where 6 is the angle between the two surfaces at x. 

On an edge V', the tangent vector ■(/', the unit normal rij 
(short for n^r. ) and the co-normal b^- form an orthonormal 
basis for each i. The circular sign convention implies the 
following: 

Lemma 2 On a Plateau border ip the first equilibrium 
equation ^ is equivalent to 

3 

^7,n, =0. (8) 

In turn, @j implies 

E-=i7,fcn(^,^,)-0, (9) 
E'=i7,fc3(^,.^,) = 0. (10) 

Moreover the pair 0), fJfl)) is equivalent to 0) at all points 
where ip is curved. 



^ The surfaces are nowhere tangent to each other. 



3.2 Proofs 

3.2.1 Proof of Lemma U 

In the plane perpendicular to if at x, the two orthonor- 
mal bases (nF,n_F A (f) and (n5,ns A (p) are rotated by 
6 with respect to each other. Equation ((JJ) is the coor- 
dinate change formula applied to the curvature vector k 
(Appendix ifij. □ 

3.2.2 Proof of Lemma Q 

The vector product of with ip gives (jH))- 

One gets (jSJ and H1Q|I by taking the scalar product 
of with, respectively, the curvature vector h = ip and 
"0 A k and applying the definitions of normal and geodesic 
curvatures (|42|l . When k / 0, the reciprocal is also guar- 
anteed because (k, A k) constitutes an orthogonal basis 
of the plane perpendicular to ip where the vectors in ^ 
are always confined. □ 

3.2.3 Proof of equation Q 

Because of normal incidence, the tangent (pj (0) coincides 
with the co-normal hj := h{ip, Tj) at p. Thus Q is just a 
rewriting of (j^J- □ 

3.2.4 Proof of equation © 

By normal incidence again, ^ and ipj form an orthonormal 
basis in the tangent plane TpTj. Therefore, as trace of the 
curvature operator H41|l . 

The linear combination with coefficients 7^ yields 
3 

= ^3 ^" ' -^J ) + X! ^" '^'^ ' ^3 ) • 
i=i J 3 

The left hand side is zero by Laplace (PJ. The last term 
on the right vanishes by Normal incidence, (TJ with 
^1 = 7r/2, gives kg{(pj,S) — —kn{(pj,J-j) and concludes. 
□ 



4 Conformal invariance 

Conformal invariance is a remarkable property of the mod- 
els of equilibrated foams. Of course, the exact symmetry 
group depends on the dimension and on the model, that 
is, on the set of equations fulfilled at equilibrium. The 
equations may be divided into two sets: i) the film, or cell 
boundaries, equation ^ (Laplace- Young); ii) Plateau's 
laws (j21El): or their 2D counterparts lOEllj for the Plateau 
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borders, which are the co-dimension 2 skeleton of the cel- 
lular complex. 

In standard foams, without any externally applied field, 
the films are constant mean curvature surfaces (CMCS). 
In general, this property is not preserved by conformal 
mappings of the embedding space or region. 

When external fields, such as gravity or electromag- 
netic forces, act on the foam, the films may not be CMCS 
anymore. This is also the case when the embedding space 
is not flat, as for a foam on a curved wall. Indeed, the pres- 
sure difference is the sum of two curvatures: the geodesic 
(longitudinal) one and the normal one, which may both 
vary along the curve even if the sum H is constant. This 
is the more general situation considered here. If only the 
second set of equations is kept (Plateau), then the sym- 
metry group is the group of conformal transformations. 
Let us treat the higher dimensional case first. 

4.1 Conformal invariance of the Plateau borders 
equilibrium. 

In Euclidean spaces of dimension larger than 2, any con- 
formal transformation is the composition of similarities 
(translations, rotations, dilatation) and inversions (Liou- 
ville Thorin theorem 13 ). 

In a dry 3 or D-dimensional foam, consider a triple 
made of a Plateau border ijj and three films jFj, j = 1, 2, 3. 
Then the line equilibrium equations 0) and |^ are pre- 
served by conformal transformations of space. 

Our proof is similar to that of By definition of 
conformality, equation |(2Jl, stating that the dihedral angles 
between the films J-j at ip are constant, is invariant by 
conformal transformations. 

Next, let us treat 101. In Liouville Thorin's decomposi- 
tion, the similarities obviously preserve and Q, so we 
only need to check invariance for inversions. Let Vj € M'^ 
be the position of a point in J-j and Hj the mean curva- 
ture at that point. Under the inversion x = the mean 

curvature transforms to Hj given by [T^ : 

H, ^ -H,r'^ - 2 ryn,. (11) 

Then for any point r of the image of ^p, thus common 
to all the J-j, the linear combination with coefficients 
yields 

3 3 3 

j=i j=i j=i 

The terms on the right hand side of H12() vanish by ^ and 
JSJ) respectively, proving the claim. □ 

Obviously, the tetrahedral figures around equilibrated 
vertices in 3D space are conserved by angle preserving 
transformations. 

In 3 and higher Z3-dimensional foams, the conformal 
group of similarities and inversions is a symmetry of the 
entire foam, including the films satisfying Laplace's equa- 
tion, only if the films are spherical caps; this is a very 
restricted subclass. 



In the 2D Euclidean plane, however, any CMCS is a 
line of constant curvature, that is, a circle. Then, indeed, 
the symmetry group of 2D foams is SL(2,C)/{1, —1}, the 
group of linear fractional transformations (or homogra- 
phies or Moebius transformations) |8I9I1U| . The two sets 
of equations, including Laplace, are preserved in this case. 
Foams in standard flat Hele Shaw chambers fulfll all these 
conditions. 



4.2 Conformal invariance of bi-dimensional foams 

Specific to 2D is the existence of a much larger set of 
conformal transformations. This larger set is a symmetry 
of foams as characterised by the Plateau set of equations 
only. The geodesic curvatures kj are not necessarily con- 
stant, the graph edges are of any smooth shape, but nodal 
equilibrium (Plateau conditions) is still assumed. 

A further extension is that the degree (or coordination) 
d of the vertices is arbitrary (but finite) , whereas standard 
2D foams have degree d = 3. 

The arguments leading to conformal invariance are de- 
tailed in |14) . For consistency, the main steps are reported 
here. The embedding plane may be identified to the com- 
plex plane: ~ C. In the plane, a foam is a planar graph. 

4.2.1 Equilibrium in 2D 

The 2D equilibrium conditions ((Sj 13 at the vertices may 
be rewritten as 

E7,T, =0, (13) 
E7,fc,=0, (14) 

where kj are the film geodesic curvatures and tj the unit 
vectors tangent to the films. The sums are over the d edges 
incident to the considered vertex. 



4.2.2 Conformal maps 

Two dimensional conformal maps are given by complex 
holomorphic functions / of the variable z — -\- ix^ in 
C ~ M^. Indeed, Cauchy-Riemann's equation, 

5.*/ = ^(ai+*a2)(/i + */2) = 0, (15) 

is equivalent to the conformality condition on the Jacobian 
matrix. The complex derivative is written /' = dzf ■ 

The scalar product in the plane is x ■ y = Re(a;*?/) = 
Re(a::y*) for any complex x,y. 

4.2.3 Paths under conformal maps 

The edges of the graphs are described by parametrised 
curves ip : [ti,t2] (or C), 1 1-> (pit). The definitions 
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of the local basis (r, n) and curvature k are recalled in 
appendix 1X1 

By a conformal /, a path ip is mapped to an image 



If : [ti,t2] 3 t (p{t) 
formation rules are 



/ 



!p{t) = f{(p{t)) e C. The trans- 



l/'l 



l/'l 



1/1 



Re 



./" ■ 



(16) 
(17) 

(18) 
(19) 



where /' = f{^{t)), etc. 



4.2.4 Conformal invariance of nodal equilibrium 

Let X — (y,i?,7) be a planar graph with vertex set V, 
edges E (given as twice differentiable curves) and line ten- 
sions 7 : i? — > M. 

Let f be a conformal map of the plane. Then the equi- 
librium conditions ]1SI are satisfied at all vertices of 
X ~ fix) if and only if they are satisfied by X . 

Indeed, any vertex x of the graph is a point common 
to all its incident edges j = 1, . . . ,d. 

The equivalence for is an immediate consequence 



of conformality as, by (fll^ . J2 Ij 
For (ini), we may use ((T^ : 



l/'(^) 



d 



1 



2^ ij i^ 'j 



Re 



f" 



(20) 

where everything is evaluated at x. By l|13|) . the real part 
on the right vanishes so that the two sums, around x and 
X, are proportional. □ 

Drenckhan et al. gave an alternative argument, based 
on the known invariance of Plateau's laws under homogra- 
phies and Taylor expansion to second order. 



4.2.5 Generalisation 

The embedding plane C may be replaced by analytic (Rie- 
mann) surfaces. Invariance and the first proof still hold 
true in these more general cases where X is a foam in a 
first surface M, X = f{X) is the image of X by a confor- 
mal map to a second surface: f : M ^ M. 



5 Mapping the edges: foams between non 
parallel plates 

In the previous section, we have seen that the equilibrium 
conditions at the Plateau borders were invariant under 
conformal maps. Here, on a specific experimental set-up, 
we consider the films more closely. 




h(x+dx) 
h(x) 



X x+dx 

Fig. 2. A film joining two non parallel plates. 



Among all the planar conformal maps, only the homo- 
graphies, or Moebius transforms, map circles to circles. In 
the traditional setting of the Hele-Shaw cell, normal inci- 
dence to the plates implies that the film curvature in the 
vertical direction (normal to the plates) is zero. Whence 
the mean curvature H reduces to half the geodesic cur- 
vature of the trace of the films on the plates, which is a 
line. By Laplace, constant pressure in each of the bubbles 
implies that the geodesic curvature is constant along any 
bubble edge, therefore the edge is a circular arc. 

Under mappings such as the logarithm or any other 
one not equal to a homography, the image of a circle arc 
is not a circle arc. So the resemblance of the pictures with 
the complex map may be questioned. 

When the foam is sandwiched between plates that are 
not parallel, the films must bend in the direction normal 
to the plates. This slightly modifies the pressure balance 
as we now examine. 



5.1 Thin foam between two non parallel plates 

The upper plate — more precisely its face wet by the foam — 
is described by a height function h{x^,x'^) defined in the 
plane at the bottom of the cell. Both h and its slope |Vft.| 
are supposed to be small. As a two dimensional vector, 
V/i locally indicates the direction of maximal slope. The 
slope in this direction is the local dihedral angle (which 
is small). Because the film meets both top and bottom 
faces orthogonally, it must be curved: to lowest order, its 
sectional, or normal, curvature in the vertical direction is 
fc„ ~ —njr ■ Vh/h (fig. El). Within the same degree of ap- 
proximation, we get the mean curvature H by adding the 
horizontal sectional curvature, which nearly is the geodesic 
curvature k := kg{ip,S) of the curve ip = J-'nS in the up- 
per or lower boundary S. Then, Laplace's law says 



Pi - Pa = = 7 (fc -I- fc,) ~ 7 fc 



(21) 



Thus, in equilibrium, with pressure difference constant 
along the film, the geodesic curvature differs from a con- 
stant by the last term in H21|l . representing contributions 
from curvatures in the third direction. 

Next, let us explicit the deformation from circularity. 
Suppose the foam X is conformally related to a reference 
one hy X — f{X). Any edge (p has an image p — f o ip in 
the reference X. Using n = ir, let us solve formula (|19|l 
of for the curvature vector k: 



|/'|fc-|-Re[n(ln/')'] 



(22) 



6 M. Mancini and C. Oguey: Foams in contact with solid boundaries: equilibrium conditions and conformal invariance 



In our approximations, njr ~ n. Fusing 1)21(1 . translated 
into complex notations, and ((22|l gives 

(Pi -P2)/7^ l/'|fc + Re[n9,(ln/'-21n/i)]. (23) 

This equation, valid on any edge of the foam, is a con- 
straint relating the conformal map / and the cell profile 
h] it involves the edge through the pressure drop, the cur- 
vature of the reference edge and the unit normal. 

For an equilibrated foam containing gas, the pressure 
drop on the LHS should be constant along any edge. Now, 
except for special cases, making the RHS constant seems 
to imply an intricate relation between the geometry of the 
edges, the mapping / and the profile h. 

Interesting cases are when the mapping / leaves some 
freedom to the choice of the foam. For example, global mo- 
tion by translation and rotation, eventually with appro- 
priate deformation. These criteria are met when Pi — P2 
and k are both identically zero. Let us treat some of these 
cases. 



5.2 Constant pressure foams 

For specificity, we assume that the reference foam is in a 
standard parallel Hele Shaw cell. Then the curvature fc, 
on the right, is proportional to the pressure difference in 
the reference: Pi — P2 = jk. 

When, moreover, the pressure is constant in both the 
curved and reference foams, equation 1(23(1 reduces to 

Re[nd,{\nf' -2\nh)] ~ 0. 

If the bubbles are small compared to the length scale over 
which / and h vary, then the 2D vector dz (In /' — 2 In h) 
must vanish. As In /' is analytic and In h real, a little rea- 
soning leads to the condition \nh/ho = Re [In/'], or 

h/ho ^ l/'l (24) 

where ho is positive a constant. 



Examples 

Following Drenckhan et al. [7! , we will most often take the 
hexagonal foam as reference (fig. EJ- 

5.2.1 Log map 

The logarithmic map 

FiS) = ^\n{a*~z), (25) 

where a is a complex number or 2D vector, maps the out- 
side of the unit disc to the upper half plane when a is 
purely imaginary. This is one of the examples experimen- 
tally demonstrated in 7 . This mapping is also involved 




Fig. 3. The reference regular hexagonal foam. 




Fig. 4. Image of the regular hexagonal foam by the log map. 

in gravity arches observed in ferro-fluid suspensions un- 
der gravity [15I16[ and in ferro-fluid bi-dimensional foams 
(17j . It was observed that 'In' is the only conformal map 
translationally invariant [15llf)ll8j . 

Applying to the inverse mapping /(z) = F^^{z) ~ 
exp(a*z), the solution for the profile is 

h{z) oc l/'l = exp(Re[a*z]) — exp(a • z). (26) 

This exponential profile matches the experimental one, at 
least within the observable precision. 

With formula ((22(1 . we can explicit the edge curvature: 

k ~ |a<^|fc + a ■ n. (27) 

In ((27(1 . the first term on the right is the effect of the 
local scale change by |/'| (A; = in a constant pressure 
reference, as assumed here) ; according to the second term, 
the edges nearly parallel to a are almost not deformed, up 
to similarities, whereas those perpendicular to a are bent. 
See fig. ^ where a (x i is vertical. 

5.2.2 Power laws 

Upper plates with circular symmetry are described by 
height functions depending on distance only: h — h(\z\). 
For the conformal map, the candidate functions are of the 
form /(z) oc z^/a. An exponent a = m/6, m integer, 
will induce a disinclination of the hexagonal crystal. The 
isobar condition 1(24(1 then implies 

h{\z\)^\f'iz)\^\zr\ (28) 



M. Mancini and C. Oguey: Foams in contact with solid boundaries: equilibrium conditions and conformal invariance 7 




Fig. 5. This 18-fold symmetric pattern maps to the hexagonal 
foam of fig-Olby f{z) oc z"^. 

For example, a spherical profile corresponds to a — 1 — 2, 
that is, to a figure with 18-fold symmetry if the reference 
foam is the regular honeycomb (fig. EJ. This result, de- 
rived from constant pressure, differs from the prediction 
of constant volume (|29|l . which is 2 (a — f) = 2, implying 
12-fold symmetry. 

5.3 Discussion, remarks 

Equation (|23|l insures compatibility with Laplace law for 
films in a thin curved interstice. It is interesting to com- 
pare it with other conditions. 



5.3.1 Volume conservation versus constant pressure 

In [7], the set-up was similar to the conventional Hele- 
Shaw cell, except that the upper plate was not parallel to 
the bottom one. It was either tilted or moderately curved. 
The major constraint, on the mapping, was volume con- 
servation in each bubble. In the long wavelength approxi- 
mation, equal volume of the reference and the transformed 
bubbles amounts to the following condition: = 
h\f'\~'^ = const., or 

h/ho = l/f (29) 

for some constant ho- In [7|, the constant volume predic- 
tions do not always match the experiments^. 

Compared to H29|l . equation l|24|l appears as a con- 
stant perimeter condition: at constant film thickness, the 
amount of liquid per bubble would be the same in the 
reference and in the transformed froth. This is just an 

^ The clearest evidence is with a spherical upper plate, where 
constant volume predicts 12-fold symmetry (a = 2) whereas 
the sample has 9-fold symmetry. Possible explanations of this 
discrepancy may rely in boundary effects, both near the ori- 
gin where the liquid occupies a significant region and at the 
outer boundary where the absence of the outer part may in- 
duce relaxations, and departure from the ideally infinite model. 
The experimental protocol followed from the reference to the 
curved foam may be of importance too. 



interpretation because (|24|l was derived from a constant 
pressure hypothesis and equilibrium. 

As already noticed, the predictions deduced from equa- 
tions H29(l and (|24|l are different, except for the trivial case 
where |/'| = 0, corresponding to changing the chamber 
thickness while keeping the two plates parallel, and where 
(|2ni and are compatible. 

Let us see the consequences of assuming both constant 
volume and local film equilibrium. Inserting H29|l into the 
expression (|2f |l for the vertical curvature gives 

ifc^ ~ -Re[ndAiih] = -Re [ndAr^lf'f] 

= ~Re[ndA-af'] = -k (30) 

by H22I) . Therefore the pressure difference (|2f |) becomes 

iPi~P2)/l^k + k,^-k, (3f) 

a result opposite to what a direct readout from the 2D 
foam and equations would predict. In figure^ the pressure 
decreases as one moves up. 

5.3.2 Short versus long range 

fn many cases, a single bubble or liquid bridge between 
non parallel plates is unstable jJHI- So, for foams, stability 
is a collective effect; the inner bubbles repel the outer ones. 

In the long wave length limit, we can clarify the re- 
spective status of the previous equations. On one hand, 
with global symmetries such as rotation or unidirectional 
translation, bubble volume preservation, equation (|29() . es- 
sentially forces the positional ordering of the cells with 
respect to each other, implying a long range order. On 
the other hand, equation (|23|l . or its constant pressure 
version (|24|l . represents the equilibrium condition for the 
films; this is a shorter range constraint because the films 
are small, of the order of the bubble size, and bounded 
by the Plateau borders. How this effect propagates across 
film junctions is not clear yet. Neither is it obvious, a pri- 
ori, that the two conformal maps, in (|29|l and (|23|l . should 
be the same, despite the fact they were both called /. 

5.3.3 Analogy with foams in gravity 

When the film is subjected to a force field iig per unit 
surface (length in 2D), the force balance in the normal 
direction is [50] 

iPi-P2)-2jH + ^ig-n = 0. (32) 

The force balance in the tangent direction implies that the 
surface tension slightly varies with position on the film: 
7 = 7(r)- When, moreover, the fluids on both sides of the 
film are liquids of mass density pi, P2, the local hydrostatic 
pressure difference is (Pi — P2) — (pi — P2) g • r-|- const. . 
In particular, if the liquids are the same, as in ferro-fluid 
foams !T7|, or if the cells contain gas of negligible weight, 
then (Pi — P2) is constant along the film, with a good 
accuracy. 

Comparing H23|l and H32|l shows that a non parallel 
profile affects the foam as an effective force field Vz (In /' — 
21n/i). 
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6 Conclusion 

When a foam is in contact with a smooth sohd surface, the 
2D foam formed by the contacts satisfies the 2D Plateau 
equihbrium conditions in the case of normal incidence. 
The case of oblique incidence is more involved and no 
simple separation appears to occur of the 2D from the 3D 
problem (Appendix [nj ■ 

The equilibrium equations for the co-dimension 2 ob- 
jects (nodes in 2D, Plateau border lines in 3D) are confor- 
mal invariant. This has been proved in jl4j and here for 
foams in Euclidean (flat) spaces. 

On the other hand, the equilibrium equation for the co- 
dimension 1 films or interfaces, namely Laplace- Young's 
law, is in general not preserved by conformal mappings. 
For flat 2D foams, the symmetry group of equilibrium 
patterns is the group of homographies, generated by Eu- 
clidean similarities and inversion. In 3D, where the group 
of conformal transformations reduces a priori to similari- 
ties and inversions, Laplace equation is preserved only in 
the special cases where the cell boundaries are spherical 
caps. 

As an application of the equations at the boundary, we 
have analysed the films equilibrium conditions for foams 
enclosed between two non parallel plates. In the limit 
where the thickness, slope and typical wave vectors of the 
profile are small, we deduced an approximate map-profile 
equation which was then solved in the case of constant 
pressure and illustrated in conditions comparable to avail- 
able experiments. 

In most of the experimental demonstrations of confor- 
mality (non parallel plates, spherical cell a la Hele-Shaw 
or ferro-fluids under gravity) , the length scale of the map- 
ping is large compared to cell diameter. As a small scale 
condition, our equation completes longer range ones such 
as the constant volume condition. Moreover, the effect of 
slightly non parallel plates can be interpreted as an effec- 
tive force field on the 2D foam. 

Examining the projected foam only can be misleading. 
For example, on the bottom figure 01 which is the model 
for a foam confined in an (approximately) exponential cell 
as on figure 3(b) of [7], the pressure seems to increase as 
one moves up (the pressure is higher on the convex side 
of a curved edge). But this is wrong. In these conditions, 
solving the equilibrium problem in 3D, taking into ac- 
count the vertical curvature implied by normal incidence 
between non parallel plates, leads to the same image (fig. 
0)) even at constant pressure; then, at each edge, the cur- 
vature in the vertical direction is exactly opposite to the 
horizontal one seen in projection, so that AP/j = if ~ 0. 
The same computation at constant volume even predicts 
that the pressure decreases as one moves up in figure 0] 

Our derivation of the map-profile relation relies on an 
approximation where the vertical curvature is retained 
only to the lowest order, constant. Further investigations 
are still needed to control this approximation or to get 
exact results on foams between non parallel plates. 

Experiments at constant pressure, rather than con- 
stant volume, would be interesting. So would be pressure 
measurements in the experiments at constant volumes. 



In the language of defects and Volterra processes, a 
mapping of the type z i— > z", a rational, represents a 
disinclination. Now the nearly conformal foams with ro- 
tational symmetry also have dislocations, which would be 
worth taking into account. 
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Appendixes 
A Planar curves 

For a parametrised curve (p : [ii,t2] ^ (or C), 
t 1-^ </?(*), arc length is ds = \d(f \ = \ip\dt and the local 
Frenet basis (r, n) is composed of t = and n = ir; i is 



the imaginary unit. The scalar curvature is 



and the curvature vector 

dr (p 
kn = — — = — 



(^) T 



ds \ip\^ 



(33) 



(34) 



B Curvatures 

We briefly recall the basics on curvatures and curves in 
surfaces 

Up to second order around a point p the surface S is 
the graph of a quadratic form g(a, a) defined in the tan- 
gent plane TpS with ordinate axis the normal line IRns(p); 
ns is the unit normal to S. The corresponding symmetric 
bi-linear form defines the curvature (linear) operator Q by 
q{a, b) = a • Qh for all a, b in TpS. 

The principal curvatures Ai , A2 and directions Ui , U2 
are the eigenvalues and orthonormal eigenvectors of Q. 
The mean and Gauss curvatures are, respectively, H — 
itrQ = (Ai + A2)/2 and G = detQ = A1A2. 

In a local parametrisation chart {z^,z'^) ^ r(z^,z^), 
the matrices of the first and second fundamental forms are 
usually taken as 

I : r^dz* • rjdz^ = gijdz^dz^ g.^ = • rj; (35) 
II : q{ridz^, Vjdz^) — qijdz^dz^ qij = r^j • ng, (36) 

where = dv/dz^, i,j = 1, 2 and summation is implicitly 
assumed over repeated indices. If denotes the matrix 

defined by Qrt = rjQ^, then = • Qrj = gtkQj or, as 
2x2 matrices, 



(37) 

The identity • a = for all a in TpS implies ns^i — 
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The sectional curvature of S in the direction u G TpS 

is 

fc„(u,5) = II(u,u) = 9(u,u) (38) 

Completing u into an orthonormal basis (u, v) and calling 
a the angle between Ui and u, the change of coordinates 
implies (Euler's Theorem) 

fc„(u, 5*) = u ■ Qu — Ai cos^ a + A2 sin^ a, (39) 

or, since Ai, A2 are the solutions of the characteristic equa- 
tion A2 - 2HX + G = 0, 

fc„(u, S)^H+{H^ - G)i/2 cos(2a). (40) 

Finally, recall the invariance of the trace: 

2i/ = A1 + A2 = u-Ou+vOv = fc„(u,S') + fc„(v,S') (41) 

for any orthonormal pair (u, v) of tangent vectors. 

When a curve tp (parametrised by arc length s) is em- 
bedded in a surface 5, its tangent vector ip{s) is, of course, 
in the tangent plane T<^(s)iS'. Its curvature vector k = if{s), 
always perpendicular to the curve, may be decomposed 
into 

k ^ kg + kn kg{ip, S) ns A if + K (</?, S) ns (42) 

where hs — A (ys is the co- normal. The component 
kg{ip,S) tangent to the surface is the geodesic curvature 
of if in S; the other is the normal curvature kn{(p, S) which 
coincides with the sectional curvature of S in the direction 
(p: kn{ip,S). 

C Bubble-glass contacts with difFerent 
surface tensions 

When the surface tension along the solid surface differs 
from cell to cell, the incidence is not normal anymore. Let 
us treat this case. 

C.l Equilibrium equation for non normal incidence 

Consider, again, an equilibrated triple {ip, ^1,^2, ^3) Each 
of the contact angles 6j = Z{nj,ns),j = 1,2,3, satisfies 
Young's law the force balance projected tangentially 
to the surface. The component along the normal axis ng 
involves the normal reaction of the glass, which is an ad- 
ditional unknown. 

If ipj is the contact curve J-j Cl S, ipj is normal to both 
rij and ns (fig-Ellj so that 

ns ■ Uj = cos9j, ns A rij = ipj sm{9j). (43) 

The following equations follow from taking the dot and 
cross products of ns with Plateau's law ||SJ: 

E'=i7.-cos(f?,)=0, (44) 
E-=i7,V^,sin(0,) = O. (45) 




Fig. 6. Relevant vectors in ipj± (left) and in TpJ-j (right). 

Regarding the curvatures, we shall prove the following: 
at the triple point p = 'ipnS, 

y^7j [kn{(pj,S)cosej - kg{(pj,S)smej] 
3 

^J2^AH]-G,y/'cos{2a,). (46) 
j 

It involves the angle aj = Z{ui{Tj),^pj) between the 
curve (fj and the first principal direction of J-j . 

Equations (03, are the analogous of (jSJ, In 
the present case where the films are not perpendicular to 
the solid surface, the equations in the boundary do not 
decouple from the full problem of the 3D foam. 



C.2 Derivation of the equilibrium equation (|46|l 

Let (ui, U2) = {ui{Tj), U2{J-j)) be the principal directions 
of Tj. We use the following angles: aj = Z{ui,ipj), f3j = 
Z(ui,^), e.j = Z(n5,nj); 

I ■ ;\ I \ ( cos a," cos/3,\ , 
(^,.^) = (urU2)(^^.^^j^.^gj. (47) 

As both ipj and are curves in Tj , equation H4U|) reads 

fc„ ((^j ,Tj)= Hj + {H] -Gjfl'^ cos(2aj- ) , (48) 
fc„(^, Tj) = H, + (i/| - G,)i/2 cos(2/3,). (49) 

Putting the coefficients 7j and summing 1)481 149|) over the 
triple around ij: gives 

^ 7, K (^, , -^j ) = E ^1 - ) ) , (50) 

'j j 

iMi^^ ^j) = E ^^ (^1 - c°«(2/3,) = 0. (51) 

3 3 

The second vanishes by 1^. In the first, ^ allows to ex- 
press the LHS in terms of curvatures in S: 

J2jl3 [kni'fj^S) cosOj - kg{ipj,S) sinOj] 

= E,7,(^'-G,)'/'cos(2a,). (52) 

This is (gni). □ 
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C.3 Remark 

The equations for the case of normal incidence follow from 
the general ones. Indeed, \9j\ = 7r/2 = \fjj — aj\ imphes 
cos(2/3j) = — cos(2aj); then, in the cosOj term dis- 
appears from the LHS and the RHS is zero hy (j5T|l . so 
that I^Hjl reduces to 
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